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In this paper, we present the higher order spectra of a scalar field produced through the higher 
derivative interactions in the initially anisotropic universe. Although we ignore the backreaction 
of the scalar field on the geometry, our analysis should have much overlap with the quantum fluc- 
tuations of the inflaton field in the anisotropic universe. We also include the planar modes whose 
O^l ■ momenta are along the plane which is perpendicular to the primordial preferred direction, for which 

effects of the initial anisotropy are not suppressed. The presence of a negative frequency mode 
produces features distinguishable from the case of the de Sitter inflation. We also show that richer 
features appear in the trispectra due to the primordial anisotropy. 

>v 

j3 ' PACS numbers: 98.80.Cq 

1 I. INTRODUCTION 

(N : 

Recent measurements by the WMAP satellite [l|, 0] have suggested that the observed map of CMB anisotropy 
1 is almost consistent with the Gaussian and statistically isotropic primordial fluctuations from inflation. However, 
r \ [ the issues on several anomalies in the recent data of the large-angle CMB temperature map have been controversial 
■ (see, e.g., [3|-l6|). These observational facts have motivated us to investigate the possibilities of the preinflationary 
' anisotropy [7H12|. The cosmic no-hair theorem ensures that, in the presence of a positive cosmological constant, an 
Oh i initially anisotropic universe exponentially approaches the de Sitter spacetime under the strong or dominant energy 
q condition (l3j . Therefore, the initial universe may be highly anisotropic, and its consequences may be inherited in 
?-h , the cosmic observables today, such as CMB anisotropy, if the number of the e-foldings for inflation takes a minimal 
r/2 ■ value 043- ^ u an expanding Kasner universe, one of two polarizations is coupled with the scalar mode, while the 
| other is decoupled, leading to an asymmetry between them. There could be a difference between the amplitudes of 
gravitational waves due to the primordial anisotropy fill ] . 
t—{ ' In this paper, first, we will investigate the bispectra of a scalar field in the initially anisotropic universe approaching 
J> . the de Sitter spacetime (see, e.g., [T3, [lj| for reviews in the case of the standard inflation). Amoung two planar 
branches of the expanding Kasner spacetime, the adiabatic vacuum can be well defined in the branch where the initial 
spacetime structure is the product of the two-dimensional Milne universe and two-dimensional Euclidean space. Thus, 
the inflaton fluctuations arc quantized in an adiabatic vacuum state different from the Bunch-Davies vacuum [l6j], 
which could give a signature distinguishable from the primordial anisotropy. Another important effect from the 
*n ■ primordial anisotropy is the mixing of the scalar and tensor metric perturbations in terms of the three-dimensional 
symmetry [7rlTo|. In particular, in the nonlinear perturbations such effects may be more significant. The purpose of 
this paper is to clarify the effects of the change of the initial vacuum state onto the higher order spectra in a single 
scalar field theory [Hj, and hence we will not discuss the coupling effects in the metric perturbations. Although we 
ignore the backreaction onto the geometry, our analysis should have much overlap with the cosmological perturbations. 
We can expect that some features appear in the limit ki + &2 ~ ^3; which is expected to be sensitive to the mixture 
of a negative frequency mode (l7j . The recent work Ref. [i~2| has confirmed part of the above expectations. The 
work focused on the nonplanar high-momentum modes whose momentum vectors are not alon g th e plane where the 
effects of the negative frequency mode are suppressed due to the adiabaticity parameter [Tol. We expect that 
the nonnegligible contributions of the primordial anisotropy to the bispectrum would appear more significantly in the 
planar modes whose momenta are along the plane. Thus, our analysis includes them. We will also investigate the 
trispectra in which the contributions from the primordial anisotropy may become important even for the nonplanar 
modes. 

The construction of this paper is as follows: In Sec. 2, we review the framework to analyze the higher order spectra 
in the anisotropic universe. In Sec. 3, we introduce the Kasner-de Sitter spacetime which exhibits an isotropization 
due to the presence of a positive cosmological constant. We also review the quantization of a free scalar field in 
this spacetime. In Sec. 4, we compute the bispectra in the presence of the higher order time-derivative interactions. 
Similarly, we investigate the trispectra in Sec. 5. In Sec. 6, we shall close this paper after giving a brief summary. 
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II. SCALAR FIELD INTERACTIONS IN THE ANISOTROPIC UNIVERSE 

A. Canonical formulation 

In this section, we present the framework to analyze the higher order spectra of a scalar field in the anisotropic 
universe, based on the so-called in-in formalism (l8j | . 

We consider a scalar field in a curved spacetime background 

d*xV=gC = -[d i xV=g[±(d<l>) 2 + v], (1) 

where V represents the interaction terms which are of higher order in terms of the scalar field and its derivatives. The 
form of V will be specified in Sec. III. C. As the background spacetime, we consider an anisotropic universe with the 
planar symmetry 

ds 2 = -dr 2 + a{r) 2 dx 2 + b{r) 2 (dy 2 + dz 2 ) = -e 6a{t) dt 2 + a(t) 2 dx 2 + b(t) 2 (dy 2 + dz 2 ) , (2) 

where e a = (ab 2 )i is the avegaged scale factor and dt = dr/e 3a . We call the (y, z)-plane the planar direction, and 
simply denote the components of any vector in this direction by the symbol "_L." We assume the universe approaches 
an isotropic one, a — > b, in the late time limit. Then the action of the scalar field Eq. (TTJ) reduces to 



~ 6a i 2 + i^L + i</> 2 , + 2V 



a 



2 



b 2 



(3) 



where wc have defined the short-hand notation as <fi := </>. t and <fi j_ := <f>- + <fi z . Wc use the coordinate t to label time 

rather than the physical time r. The conjugate momentum to the scalar field is given by -k^ = ^£ = </>. The scalar 

field and conjugate momentum satisfy the equal time commutation relations [4>(t, x 1 ), TTj,{t, y l )\ = i5^{x l — y l ) and 
[(j)(t,x l ),(j)(t,y' 1 )] = [n^t, x z ), 7r^(t, y 1 )] = 0. The Hamiltonian density is defined by 

H = ^ - C = ^ (e- 6a ^ + ^ + + 2V) . (4) 

The Hamiltonian is given by H\cf>(t), 7r^(i); t] := J dPxH. The equation of motion for the scalar field is given by 

^ = i[ff^(i),7r (t);t]^], K 4 , = i[H[<t>{t),-K <tl {t);t}^}. (5) 

The solutions at t can be expressed in terms of a similar operator at a very early time to through the unitality 
transformation 

(j>{t) = [/^(Mo^oMMo), n^(t) = U- 1 (t,to)TT i> {to)U(t,t ), (6) 
where U(t,to) obeys the differential equation 

jU{t,h) = -iH[4>{t)^{t)-t]U{t,h)- (7) 

The initial condition is given by U(to,to) = 1. Let us briefly discuss the choice of the initial time to- In cosmology, the 
quantization of a scalar field is usually performed when the corresponding mode is well inside the horizon, and it is 
reasonable to choose to — ¥ -co, when the effects of the cosmic expansion are completely negligible. But in any model 
where the background spacetime is modified from the standard one, the ordinary description of the mode functions 
becomes valid only after some critical time fl2l [l7j. In the case of an initially anisotropic universe, we will choose to 
to be the time when the universe is isotropized, and the effective theory description becomes valid. We will discuss 
this in more detail after specifying the interaction terms. 

B. Interaction picture 

In calculating U(t,to), we decompose H into the (quadratic) kinetic part Ho and the interaction part Hi that are 
of higher order in the scalar field amplitudes, 

H[<j>(t\*4{t);t] = H Q [<j>{t)^{t);t] +ff/[0(t) > 7r^(t);t] I (8) 
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where 

fl> [>(*), 7i*(t);t] := e 6a J d 3 xV. (9) 

We will calculate U as the power series in Hi. We introduce the interaction picture and define the interacting operators 
(j) 1 (t) and 7r{(£) whose dependence is determined by the quadratic part of the Hamiltonian Hq, 



^ (t) = i [ff (*), ^ (*);*], ^(*)1 , ^ (f) = i [#«, (t) ; f] , ti£ (t) 



(10) 



where the initial conditions are given by (p 1 (t ) = </>(fo) and 7r^(£o) = ^(^o)- 

In evaluating H^cf) 1 (t), Tv^(t); t] we can take the time argument of (f) 1 and 7r^ to be any value, and we can take to so 
that 



H (t) y 4> {t);t]=H [<t>(t ), ^ (t ); t] . 



(11) 



The solution to Eq. (jTUJ) can be written in terms of the unitality transoformation 7 (t) = U a 1 (t, to)4>(to)Uo(t, to) and 
ir^(t) = U ~ 1 (t,to)TT<j,(to)Uo(t,to), where Uo(t,to) is the solution to the equation 



We obtain 



which gives 



dt 



-C/ (Mo) = -iH Q [cj)(to),ir^(to);t]Uo(t,to). 



U -\t,t )U(t,t )} = -iU^^to^l^^^to^tjU^to), 



(12) 



(13) 



U(t, t ) = U (t, t Q )T(t, to), —T(t, to) = -iHi(t)F(t, t Q ), F(toM) = 1- 
We define the shorthand notation for the interaction Hamiltonian 

H T (t) := U (t, to)Hj [(/>(t ), n^t ); t] Uq% t ) = Hj [tf(t), ir^t); t] . 
The solution for F(t, to) is given by 

F{t,to) = Texp(-i / Hi(t)dtj, 
where T is the time-ordering operator. The expectation value for some operator A(t) is given by 
(A(t)) = ( fexp(i j Hi(t)dt\ A*(t) Tcxpf- j j H^dt) V 

J to Jto 



(14) 



(15) 



(16) 



(17) 



where A(t) can be any product of <f>(t)s and Tr<p(t)s. For the evaluations of the higher order spectra, it is more useful 
to rewrite Eq. (fTT|) as 



(A{t)) = J2* N dt N 



N=0 1/10 



dt 



N-l ' ' ' 



1-2 



dt! 



Hiih), \H I (t 2 ),---\H I (t N ),A I (t) 



(18) 



III. PREINFLATIONARY ANISOTROPIC UNIVERSE AND THE BEHAVIOR OF A SCALAR FIELD 



In this section, we present the background geometry and the way to quantize a scalar field in this background. 
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A. Planar Kasner-de Sitter spacetime 

We consider the Einstein gravity with a positive cosmological constant A > 

1 



S 9 = 1T2 I d*xV=9[R-2A). (19) 



2k 2 

There is the planar Kasner-de Sitter solution with a two-dimensional Euclidean symmetry whose metric is given by 
Eq. © with 

a = sinh*(3J? r) tanlJ (|#or) , b = sinh3(3i? r) coth* (|#ot) , (20) 
where H = is the Hubble constant in the isotropized limit. The averaged scale factor is given by 

e a := {ab 2 )? = ( sinh(3 H t)) *. (21) 

The spacetime is completely regular everywhere even in the limit of the intial time r — > 0, in contrast to the other 
branches of the Kasner spacetime. The Hubble parameter in each direction is given by 

H - := - = . Jt° , (2 + cosh(3g T)), H b := \ = H Q tanh (?-H t) . (22) 

In this solution, at the initial times the spacetime structure becomes the product of a (Milne) patch of the two- 
dimensional Minkowski spacetime and the two-dimensional Euclidean space. Thus it is possible to define the initial 
adianatic vacuum state in this solution. At the later times for r — > oo, the universe approaches the de Sitter spacetime 
with the expansion Hq- In the other Kasner-de Sitter solutions which contain initial curvature singularities, the 
adiabatic vacuum cannot be well-defined in the asymptotic past [ToL ITl| . 

B. Scalar field in the Kasner-de Sitter spacetime 

Here as in Sec. II, we will work in the time coordinate dt = -fe- Note that these two time coordinates are related 
via the simple analytic relation 

sinh(-3if *) = . * ~ v (23) 

In the interaction picture, 

/r1 3 k r rl 3 k 

e* kx 0(i,k) = / e* kx (u k (t)a k + U l k WaL k ). (24) 

(27r)2 J (27r)2 

In our Kasner-de Sitter background, 

. j2 . 04/3(1,2 6fft , p\ 04/31,2 

(^ + ^) 2 )^) = 0, tf(t) := { ^ 6Ht p } = 1^ (1 - rim), (25) 
where k 2 = k\ + k^, r± := k±/k and 



£(t) := 1 - e 6Ht = e- 6a (Ve 6a + l- 1), (26) 
varies from one to zero as time t increases from negative infinity to zero. The late time solution is given by 



t-£r(-3-ffo*)3 



e »o 



«k(t) = A+u£>(t) + A k u£»*(t), u^(t) = -j= [( - 3H ty + -^j . (27) 



The coefficients A £ and A k are determined through the matching of the mode functions. In Rcfs. [13, EI, we have 
classified modes into the nonplanar high-momentum modes and the planar ones. 
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(1) Nonplanar, high-momentum modes 

The nonplanar high- momentum modes satisfy Hq <§; ki ~ fc_i_. For these modes, the WKB solutions for the early 
time can be directly matched to mode functions in the dc Sitter spacetime at the later times. The coefficients can be 
found as 



At 



k 



where Q(r±) := | — r\ contains the direction dependence appearing only in the negative frequency mode. A^ is of 

order e\ where e* is the adiabaticity parameter e := ^qI'* at the matching time t = £». Hence, we expect that these 
modes would not give significant effects on spectra and cosmic observables at the level of the power spectrum. 

(2) Planar modes 

We then consider the planar modes with Hq < k\ <C k±, which are along the constant x plane where there is a 
two-dimensional rotational symmetry. The main difference from the previous case is that there is a temporal violation 
of the WK B ap proximation during which effects of the anisotropy are encoded into the mode functions. According 
to Refs pU [Ti( , 

A + = I e -*fk-i*+M-y/%) A - = <C1 e iJ%+i*-W-yf%) (29 ) 

where q 1 := and $ := ^ r G) fc 1_ of.f^). Since e _7n?1 may be close to unity, \A\~^\ < \A^\. 



3H ° ' 3x23r(|)i? VV Ho . 

To discuss the late time behavior of the mode functions, it is convenient to work by choosing the conformal time 
dr) = e~ a dr = e 2a dt, in which the late time mode functions can be rewritten as 



«kft) = A+ui°\ v ) + A^*( V ), u£\ V ) = (-»?+!)• ( 30 ) 



In the conformal time coordinate, the averaged scale factor is given by e c 



1 



C. Interaction Hamiltonian 



In this subsection, we focus on the choice of V in Eq. (fTJ). We consider the high energy corrections to the kinetic 
terms. Although in general as such corrections there would be both the time- and spatial-derivative terms, for 
simplicity we focus on the time-derivative terms. The time derivatives of the scalar field in the Lagrangian density 
except for the volume factor can be expanded in terms of the amplitude 



n— 3 



(31) 



where q n [n = 3, 4, ■■ •) are dimensionless constants which are typically of order unity, and M is the energy scale 
at which the higher order corrections become important. It is also appropriate to use derivatives with respect to 
the proper time r, rather than another time coordinate, such as t. We may regard our scalar field theory as the 
low energy limit of a complete UV theory, appear just as low energy corrections. In the interaction picture, the 
leading order quadratic term gives the ordinary kinetic term, and higher order ones give interactions. Since we are 
particularly interested in the three- and four-point interactions, we focus on the terms of n = 3 and n = 4. We expect 
that the time-derivative interactions would give us information enough to see at least the difference of the anisotropic 
background from the de Sitter inflation. 
Thus wc choose V in Eq. (fl]) as 



V = e-^l 



(32) 



Using 7T0 = <fi and dt = e 3a dr, Eq. (|32]l is rewritten as 



V 



S3 



AP +qi \Jp) 



(33) 
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which agrees with the n — 3 and n — 4 terms in Eq. (|31| . From Eq. ([9]), the interaction Hamiltonian is given by 

Hj[^t),7r(t);t] = £T< 3 > + [<f>(t),ir(t);t] 



- 6a 1 <Px<b\ 



<P.t . (9,r\ 



HP[<j>(tU(t);t] 



X M^ ' 

3„ 94 U 



(34) 



Let us focus more on the reasons for considering the derivative interactions. The first reason is that in the metric 
perturbation theory (most of) the interactions are given by the interactions of derivatives [l4j|. The second reason 
is that higher derivative interactions may enhance the higher order spectra JT3, H(| because the interactions depend 
on the inverse power of the (averaged) scale factor, which is large in the past. Although we ignore the backreaction 
of the scalar field on the geometry, our analysis should contain overlap with the realistic analysis on the quantum 
fluctuations. Since (\4>,t\ 2 ) — Hq, the scalar field does not backreact significantly on the background geometry so long 
as Hq < M p . In addition, in the equation of motion of the scalar field the corrections of the higher order derivative 

terms to the kinetic term appear in order of Thus, the higher-order derivative interactions can always be treated 
as perturbations so long as M > H , and then the early time WKB solutions and late time de Sitter ones given 
by Eq. (|27|) are good approximations of the mode functions with a very high accuracy. Hence, our pcrturbative 
arguments are valid when the relation Hq < M < M p is satisfied. We choose the cut-off time for the time integral to 
be = M when the physical momentum for the corresponding mode becomes below the energy scale of the derivative 
interactions. This gives —770 = tjt , where k is the total momentum for the relevant spectra, i.e. k = k± + fc 2 + k 3 for 



the bispectra and k = ki + fc 2 - 
mode function becomes valid, 
scale satisfying 



&3 + &4 for the trispectra. This time becomes later than the time when the de Sitter 
-770 < — 77, := yj= , leading to k < In other words, our interest is in the energy 



M 
H 



< el 



(35) 



for a corresponding number of k. Thus, in the subsequent parts of this paper we will focus on these modes. 

In the next sections, we will present the bispectra and trispectra following the formulation presented in the previous 
section. The leading order contribution to the bispcctrum is given by the three-point contact interaction diagram. 
Similarly, the leading order contribution to the trispectra is given by two kinds of the connected diagrams. The 
first one is the contact interaction diagram [see, e.g., Fig. 1 (A) in pjj]], and the second one is the scalar exchange 
interaction diagram [see, e.g., Fig. 1 (B) in [19j]. We expect that loop corrections are suppressed and ignore them in 
the subsequent discussions. 



D. Power spectra 



Before closing this section, we show the power spectra including the leading order corrections due to the primordial 
anisotropy, for both the nonplanar high- momentum modes and the planar modes, which were obtained in [ToL [Til]. 
Note that, in contrast to the power spectra obtained from the standard isotropic scenario, they depend not only on 
the magnitude of the momentum vector k = |k|, but also on its direction. For the nonplanar high-momentum modes, 
the power spectrum is given by 



™-(§)V«'j(fr~(VS)- 



(36) 



where Q(r±) is defined below Eq. P5j) . Thus the corrections due to the primordial anisotropy are suppressed by e 3 . 
On the other hand, for the planar modes, the power spectrum is given by 



/ Hq\ 2 1 , . , sin(2$) \ 
P(k)=(^) (coth^^-ji-L), 



(37) 



where q\ and $ are defined below Eq. ([2U]) . Thus the corrections due to the primordial anisotropy are not suppressed. 
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IV. BISPECTRA 

The leading order contribution to the bispectra is given by the contact interaction 

Mf,kiMf,k 2 Mt,k3)>=i / dti([F} 3) (ti),^(t,k 1 )/(t,k 2 )^(t,k a )]). (38) 

J to 

The bispectrum in the t — > limit takes the form of 

(^(O > ki)^(O > k 2 )0(O,k 3 )) = (2^)- 3 ^B(k 1 ,k 2 ,k 3 ) ( 5(k 1 +k 2 +k 3 ). (39) 
After computations, we obtain 
B(ki,k 2 ,k3) 

J 4+A k3 J(fc l! * 2) -ks) + 2 perms) - (.4+ A^l{h, -k 2 , -k 3 ) + 2 perms) 
+ C.C.}, (40) 

where 

■ 1 , f 2z - e -*(pi+»-H»)'» [2i - 2(p! +p 2 + p 3 )% - *(pi + P2 +Pa) 2 % 2 ] }, 



4(fcifc 2 fc 3 ) 

X 



(Pi +P2 +P3)' 

and Z(pi,p 2 ,p 3 )* = X(—pi, — p 2 , — p 3 ). In the limit of pi +p 2 +P3 — > 0, the function X remains finite, 



(41) 



l(pi,P2,p 3 ) = y+0(Pi+P2+P3). (42) 

Since |Z(pi, j» 2 , j» 3 )| is increasing in the limit of pi + p 2 + p 3 — > 0, the bispectrum would have peaks in the limits of 
k\ — > fc 2 + fc 3 , fc 2 — > fci + fc 3 and fc 3 — > fci + /c 2 (case I). The bispectrum can also have a peak in the limits k\ — fc 2 3> fc 3 , 
&2 = &3 3> ^1) an d fc 3 = fci 3> fc 2 (case II). We discuss each case separately. 

1. Case I 

In the limit of of fci — > fc 2 + fc 3 , the leading order contribution to the bispectrum is given by 
B(ki,k 2 ,k 3 ) 

(43) 

Similar results are obtained in the limits of fc 2 — > ki + fc 3 and fc 3 — > k\ + fc 2 . For the bispectrum for the nonplanar 
high-momentum modes, 

Q(r±,i) Hi 



2 k\ 



- ^)*« - A U2 Y( A t ^ 3 y(A u A+ 2 A+ 3 A+ A k2 A k j} 

(44) 

Thus 



^p 3 
4(fc^ 2 A : 3) , 



S(k!, k a , kg) -> "V" r Q(r x ,i). (45) 
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The ratio of Eq. (1451) to the case of the de Sitter inflation, Eq. (|A1[) . is given by 

M 3 fri H 2 

Q(r^)Hlk lV l = Q( rXil )-^^_o , (46) 
which is of order 44" e t- Thus the direction dependent bispectrum can be as large as that in the de Sitter inflation for 

_ 4 

■j^ > e* 3 . But this is never satisfied because of Eq. (|35|) . 

On the other hand, for the planar modes, using the phase factor ip := — — $ + (j — \J jj^) the bispectrum in 
the same limit is given by 

mk k , go?jg e -^^sin(2^) 

1 J ' 2 ' 3J 2(fc 1 fc 2 /fc 3 ) (i_ e -*rgi.i)(i_e- a '9i.a)(l-e-* r ft.»)' 1 J 

The ratio to the case of the de Sitter inflation is given by 

M 3 kf e -*qi,i sin(2V>i) 



H 3 k 3 (l- e -2^l,l)(l- e -2 7 r<?l,2)(l_ e -27r gi ,3)^ 

~~<>^. 1 ' I 'Ti n n Vvi~.nn^] in nAnni ttti-J-1-i T7 1 y-i /l Q ft f-f-i-v 1 



(48) 



which is greater than unity since jj- ^> 1. This bound is consistent with Eq. ([55)) for f < < e^ 1 . Therefore, 
such a bispectrum is more important than for the case of the nonplanar high-momentum modes. This bound can be 
consistent with Eq. (|35|) . 

Similar results can be obtained for the limits of fc 2 — > fc 3 + ki, and fc 3 — > k\ + fc 2 . 



2. Case II 



In the limit of k\ = fc 2 2> fc 3 , the leading order contribution to the bispectrum from the nonplanar high- momentum 
modes is given by 

S(kl ' k2 ' k3) 

- iff (Q(r±*) + Q(r±*))- ( 49 ) 

The ratio of Eq. (|49| to the case of the de Sitter inflation, Eq. (|Al|) , is given by 

/ \ M 3 k^ H 2 / \ 

fe^(Q(rx A ) + Q(rx, a )) = ^^(Q(r±,i) + OK a )), (50) 



which is of order ttte*- Thus the direction dependent bispectrum can be as large as that in the de Sitter inflation 

_ 4 

for jf > u 3 . But this is never satisfied because of Eq. (f3"5|). In addition, assuming that k = k\ + fc 2 + k$ is fixed, 
the ratio of the bispectrum of case II, Eq. (|49[) . to that of case I, Eq. (|45|) (with the permutation to fc 3 — > k\ + fc 2 ) is 
given by Thus the bispectrum in case II is enhanced to that in case I by the ratio -r- ^> f . 
The leading order contribution from the planar modes is given by 

Ic It ^ ~ H ^ e ~* qi ' 1 sin ^ + 6 ^ qi ' 2 Sin(2 ^ 2) 
[ u 2 ' 3> ~ 2{k 2 k 3 )(l~e- 2 ^){l-e- 2 ^)(l-e- 2 ^)' [ ' 

The ratio to the case of the de Sitter inflation is given by 

M 3 k 3 e ~ nqi - 1 sin(2-0i) + e^ 1 - 2 sin(2^ 2 ) 



H 3 k 3 (l- e -2^i,i)(i_ e -2 7 r< ?1 , 2 )( 1 _ e -2^ 1 ,3)' ( 52 ) 

which is greater than unity since jj- ^> 1. This bound is consistent with Eq. (f3"5j) for 1 < < e" 1 . The ratio of Eq. 
P7)l with permutation to fc 3 — ► k\ + fc 2 to Eq. ([ST]) is given by 

k 3 e -^n,a sin(2^ 3 ) 



k e- 7r « 1 . 1 sin(2-0i) + e"^ 1 - 2 sin(2V> 2 ) ' 



(53) 



which becomes smaller than unity as for the nonplanar high-momentum modes. 
A similar result can be obtained for the cases of fc 2 = k 3 3> fci and fc 3 = ki 3> /c 2 . 
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3. Features of the bispectra 

We summarize the features of the bispectra in our background spacetime and compare them with those in the de 
Sitter inflation. In our background the bispectrum in the limit of case I (for example, for ki — > k 2 + k 3 ) is greater 
than that in the limit of case II (for example, for k\ = k 2 ^S> fc 3 ), by a factor of where k is the typical magnitude 
of the momentum vectors. 

We then discuss the direction dependences of the bispectra. Though we focus on the nonplanar high-momentum 
modes, the essence is the same also for the planar modes. In case I, for example, for k\ — > k 2 + fc 3 , the bispectrum 
is proportional to Q(r^ i) and hence depends on the direction of the momentum vector ki. In case II, for example, 
for k\ = k% 3> fc 3 , it is proportional to Q(rj_i) + Q{f_\_,2)i and hence depends on the directions of the momenta ki 
and k2, which have shorter wavelengths than k3. Thus the bispectra are direction dependent, and the dependence is 
significantly different in various limiting cases. 

We finally compare the amplitude of the bispectra with those in the de Sitter inflation. For the planar modes, the 
bispectra in both the limits of cases I and II can be greater than those in the de Sitter spacetime. On the other hand, 
for the nonplanar high-momentum modes, the bispectrum in both cases cannot exceed those in the de Sitter inflation. 
Thus, combined with the direction dependence, the bispectra for the planar modes in all the limiting cases would be 
particularly interesting to distinguish models. 

In the next section, we will investigate the trispectra. In contrast to the bispectra, even for the nonplanar high- 
momentum modes, the trispectra could be more important than those in the de Sitter inflation. 

V. TRISPECTRA 

The leading order contributions to the trispectra are given by the contact interaction as welll as the scalar exchange 
interaction. 

A. Contribution from the contact interaction 

The leading order contribution of the four-point interaction is given by 

^(tMmtM)4>{tM)4>{tM)) =ij dii([ J tfi 4) (ti),0 7 (t,k 1 )0 / (t,k 2 )0 / (t,k 3 )0 7 (t,k4)]). (54) 

The trispectrum in the t — > limit takes the following form: 

(0(0, kiMO, k 2 )0(O, k 3 )0(O, k 4 )) = (2ny 3 ^ I T c (k 1 , k 2 , k 3 , k4)5(k! + k 2 + k 3 + k 4 ). (55) 
After computations, we obtain 

r c (k 1 ,k 2 ,k 3 .k 1 ) = ajSL^ - AiS (At. - A i% )*« - A £ ,)*(A£ - A^y 



I 



{kik2k%kn 

KAt A t A t^^ k ^ fc 3, fc 4 ) - A^A^A k3 A ki J(-k u -k 2 , ~k 3 , -k A ) 
A t A t A t A ^( k ^ k ^ fc 3, -ki) + 3 perms) 
A kA : k2 A k 3 A k i ^( k ^ k 2,-h,-ki) + 5 perms) 
KA^KA^^' ~ fc 2, -k 3 , -k^ + 3 perms) } 
C.C.I , (56) 
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where we have defined 

J(Pi,P2,P3,Pa) 

24i - (pi +p 2 + P3 + Pi)r]o 



4(pi +P2 +P3 +P-l) 5 

x (24 - i(pi +p 2 +P3+Pi)m{ - 12 - (pi +P2+P3 +Pi)VO 

X (4 i -(pi+p 2 +p 3 +P4)??o)))]}. (57) 

In the limit of p\ + P2 + P3 + p± — > 0, the function J remains finite as 

J(pi,P2,P3,Pi) = -^ril + 0{pi +P2+P3+Pi)- (58) 

Thus the trispectrum from the contact interaction depends on the size of four vectors k t (i = 1,2,3,4), and not on 
the size of two independent relative vectors ki2 = |ki + k2 1 and ku = |ki + ki|, since there is no internal line in the 
diagram. Since \J^(pi,P2,P3,Pa)\ is increasing in the limit of p\ +P2 +P3 +Pa — > 0, the trispectrum has the maximum 
amplitude in the limits of k\ — > fc 2 + k 3 + &4, fc 2 — > k\ + k 3 + ki, fc 3 — > fcj + fc 2 + k&, and £4 — > fci + fc 2 + &3 (case 
I), as well as ki + fc 2 —> k3 + k±, and k\ + k± —> fc 2 + £3 (case II). We also discuss the limits of ki = fc 2 ^> ki2 and 
&3 = &4 3> ^12, an d similarly fci = fc4 ^> fci4 and = fc3 ^ fci4 (case III). 

Case I 

In the limit of fci — > fc 2 + A: 3 + ki, the leading order contribution to the trispectrum is given by 

T c (k llk2 ,k 3 ,k4) -> - 5{ ^^ {{K - 

x (KAiAA + A iA* 2 K 3 K.)\ ( 59 ) 

For the trispectrum for the nonplanar high-momentum modes, we obtain 
Q(r±,i) ^ 



(60) 



Hence 



r c (k 1 ,k 2 ,k 3 ,k 4 )^- T7T773 f 1 ^-Tg (r . ±il ). (6i) 



3go 10 ^ 5 

10(fc3fc 2 fc 3 fc 4 ) 

The ratio of Eq. (I5T1) to the case of the de Sitter inflation, Eq. (|A2[) . is given by 



klH^ 5 Q(r ±il ) = ^^QK,!), (62) 

which is of order ^ e*. Thus the direction dependent trispectrum can be as large as that in the de Sitter inflation 

for -gjj- > e* 5 . With the condition for the modes Eq. (|35|) . the energy scale for the interaction is bounded as 
-i < m_ < -1 

On the other hand, for the planar modes, the trispectrum in the same limit is given by 

T(V V VVA^ 3H ^° e -^sin(2^) 

cl X ' 2 ' 3 ' 4J 5(fcifc 2 fc 3 fc 4 ) (l-e- 2 ^.i)(l-e- 2 ^ 1 - 2 )(l-e- 2 ^.^)(l-e- 2 ^^)' 1 ' 

The ratio to the case of the de Sitter inflation is given by 

M b k\ e" 7 ^ 1 sin(2V>i) 



H% k 5 (1 - e- 2,r ?i>i)(l - e- 27r «i-2)(l - e -2^i, 3 )(i _ e -27r 9l , 4 ) ' 

which is greater than unity since jf- ^ 1. This bound is consistent with Eq. ([33)) for 1 < < e^ 1 . 
Similar results can be obtained for the cases of fc 2 — > ki + k3 + ki, k3 — > fci + fc 2 + ki, and ^4 — > ki + fc 2 



(64) 
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2. Case II 



In the limit k\ + fc 2 — > fc 3 + k±, the leading order contribution to the trispectrum is given by 

3H$r, 5 



T c (ki,k 2 ,k 3) k4) -> 



5(/ci/c 2 £; 3 fc 4 ) 



A -k 2 A k 3 A ii + A kl A i 2 Ac 4 ) } • 



(65) 



We see that this type of trispectrum can be much more important than in the previous case. For the trispectrum for 
the nonplanar high-momentum modes, we obtain 



Im « - A ki Y(Ai 2 A k2 Y(A+ 3 A k f(Al A ki 



A kl A k2 A+ 4 + A£ A+ 2 A k A k 



.)} 



0(rx,i)Q(rx,a) K 



7/2 



k\ /2 kl /2 ^kl k\ 



Q(r±, 3 )Q(r±A) Hi 



7/2 



( 1 1 

^ 1„ 2 JU 2 



Thus 



T c (ki,k 2 ,k 3 ,k 4 ) -> 



20(fcifc 2 fc 3 fc 4 ) 



Q(rx,i)Q(rx, 2 ) ^ J_ 



3 3 



, 3/2, 3/2 . 
Kq fc 4 fc 



Q(r±, 3 )Q(r±,i) ( 1 

3 3 I 1 

1,21.2 H2 



(66) 



(67) 



The ratio of Eq. (|67| to the case of the de Sitter inflation Eq. (|A2j) is given by 

-Q(rx,i)<90x >2 ) ( 1 , 1 \ , Q(r^ 3 )Q(r ± . 4 ) / 1 



(fcx + fc 2 ) 5 ^ 



^1 ">| 



k 2 



h 2 h 2 
«, 3 h. 4 







k^ 




i 


)] 


kj 





(68) 



M 5 (fci + fc 2 ) 5 ff ' r Q(r ± .i)Q(r ± . 2 ) / 1 1 N Q(r±, 3 )Q(r±,4) / 1 

EJ"5 7 1 ; 33 I 1 T 1 ) T 33 I 1 

tin. K L L2I,2 \,L2 jL2 / 1,21.2 \ 1. 2 

' n,-^ 2 3 4 3 

which is of order 775- el- Thus the direction dependent trispectrum can be as large as that in the de Sitter inflation for 
■gt > e* 5 . But this is never satisfied because of Eq. In addition, assuming that k = k\ + fc 2 + fc 3 + fc 4 is fixed, 

the ratio of the trispectrum of case II, Eq. (joTf . to that of case I, Eq. (j()Tj) . is given by f 2 j which is suppressed. 
On the other hand, for the planar modes, the trispectrum in the same limit is given by 

3H*4 



T^k^ka.kg.k*) ->■ 



5(fci/c 2 /c 3 A; 4 ) 

e -^(«i,i+«i, 2 ) sin (2^ 1 + 2^ 2 ) + e - w (9i, 3 +9i,4) s i n (2ijj 3 + 2ip 4 ) 



(1 _ e _27r<jl ' 1 )(l — e~ 27Tqi - 2 )(l — e _27r9l > 3 )(l — e~ 27Tqi < 4 ) 
The ratio to the case of the de Sitter inflation is given by 

M 5 (ki + k 2 f e - n ^< l+qi ^ sin(2^! + 2^ 2 ) + e'^ 1 -^ 1 -^ sin(2^ 3 + 2^ 4 ) 



H 5 



k 5 



(1 - e- 27r « 1 . 1 )( 1 - e- 27r 9 1 - 2 )(l - e- 27r «!.3)(l - 



,— 27rgi,4 



(69) 



(70) 



which is greater than unity since jf - 3> 1. This bound is consistent with Eq. (pt5j) for 1 < jf- < 1 . The ratio of Eq. 
([S3]) to Eq. ([55]) is given by 



e -T«i,i s in(2^i) 



-^(91,1+91,2) s in(2V>i + 2^2) + e- 7r (« 1 ^ + « 1 - 4 ) sin(2?/> 3 + 2^ 4 ) ' 



(71) 



which is of order unity, in contrast to the one for the nonplanar high-momentum modes. 
Similar results can be obtained for the case of ki + k 4 — > fc 2 + k 3 . 
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Case III 



In the squeezed limit of k\ = k 2 3> fci 2 and fc 3 = k& ^> fci 2 , the leading order contribution to the trispectrum is 
given by 



r c (ki,k 2 ,k3,k4) 



{AtK 2 + A ki < ) (< A ki + A ka <)*} 

23 

S 77 2 ?? 5 / 1 1 \ 

Vr — + — ) (0(r±,i) + Q(r± )3 ))(Q(r± l3 ) + Q(r±, 4 )). 



The ratio of Eq. (|72|) to that in the case of the de Sitter inflation, Eq. (|A2[) . is given by 



(72) 



H-f^ih + h) 5 ( 1 1 



3 3 



f — + — ) (Q(rj.,i) + Q(rx, 8 )) (Q(rx,s) + Q(rx, 4 )) 

ft 1 ftg 



m 5 H$(k 1 + k 2 r ( i i 



n UU1.5 V £.2 1,2 



f — + — ) (Q(rx,i) + 0(rx, 2 )) (QK, 3 ) + Q(rx,4)), 

\ ;„ 2 i. 2 / 



(73) 



which is of order 775- e*. Thus the direction dependent trispectrum can be as large as that in the de Sitter inflation 



for jf > e~ E . But this is never satisfied because of Eq. (|35[) . In addition, assuming that k = k\ + k 2 + fc 3 + fc 4 is 
fixed, the ratio of the trispectrum of case II, Eq. (I57|) . to that of case III, Eq. (|T2"|) . is of order 0(1), and the ratio to 
case I, (|61j) . is suppressed by e*. 

On the other hand, for the planar modes, 



T 1 c (k 1 ,k 2 ,k 3 ,k 4 ) 



1 



5(kjk$) (1 - e- ?r< Ji.i)(l - e" lr « 1 - 2 )(l - e- 7r< ? 1 > 3 )(l - e-^? 1 ^) 



x e -^iiA+n,s) S i n (2ip! + 2ip 3 ) + 6-^(91.1+91,4) S i n (2^ 1 + 2^ 4 ) 
+ e'^ qi - 2+qi - 3) sin(2V> 2 + 2^ 3 ) + e -^ 9l ' 2+91 ^ sin(2?/> 2 + 2^ 4 ) 
The ratio to the case of the de Sitter inflation is given by 

M 5 (fci + fc 3 ) 5 1 

Hi k 5 (1 - e- 7r «Li)(l - e~ 7r 9i.2)(l - e ~' r 9i.3)(l - e -^«i,4) 

x |" e - ,r (9i.i+9i.3) sin(2^i + 2^ 3 ) + e _ ' 7r(91 - 1+?1 ' 4) sin(2^i + 2^ 4 ) 
+ e-^^ 1 .^ 91 - 35 sin(2^ 2 + 2?/> 3 ) + e^ 91 ' 2 ^ 1 ' 4 ) sin(2^ 2 + 2ip 4 ) 



(74) 



(75) 



which is greater than unity since -f^ ^> 1. This bound is consistent with Eq. (|35|) for 1 < < e„ l . The ratio of Eq. 
(|69| to Eq. (dU is given by 



,—"■(91,1+91,2 



"'(91,1+91,3 



> sin(2^i + 2^2) + e _7r(9l ' 8+9l ' 4) sin(2V> 3 + 2^ 4 ) 

> sin(2^i + 2?/> 3 ) + e - 7r (« 1 . 1 +9i.'*) sin(2^i + 2t/> 4 ) 



+ e -" {qi ' 2+qi - 3) sin(2^ 2 + 2^3) + e - 7,{qi - 2+qi - 4) sin(2^ 2 + 2^ 4 ) 

which is of order unity, in contrast to the one for the nonplanar high-momentum modes. 
Similar results can be obtained for the case of k\ = fc 4 3> fci 4 and k 2 = fc 3 ^> ku. 



(76) 
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4- Features of the trispectra from the contact interaction 

We summarize the features of the trispectra from the contact interaction diagram and compare them with those in 
the de Sitter inflation. In our background the trispectrum in the limiting case I (for example, for k\ — > k 2 + fc 3 + k±) 
is greater than those in the limiting case II ( for example, for ki + k 2 — > k 3 + k±) and case III (for example, for 

k\ = k 2 ^ kyi and ^3 = ^4^ k\ 2 \ by a factor of (jfa) 2 , where k ^> H$ is the typical size of the momentum vectors. 

We then summarize the direction dependence of the trispectra. Though we focus on the nonplanar high-momentum 
modes, the essence is the same for the planar modes. In the limiting case I, for example, for k\ — > k 2 + k% + k<±, the 
trispectrum is proportional to Q(rj_ i) and hence depends on the direction of the momentum vector ki. In the limiting 
case II, for example, for fci + fc 2 — > fc 3 + fc 4 , it is proportional to a combination of Q(r± i)Q(r± 2 ) and Q( r ±,3)Q( r ±,i), 
which depends on the directions of all the momenta. In the limiting case III, for example, for k\ — k 2 3> k\ 2 and 
&3 = k± 3> fci2, the trispectrum is proportional to (Q(r± i i) + Q(r± 2 ))(Q(r± 3) + Q(rx,4)), which depends on the 
directions of all the external momenta. Thus the trispectra are direction dependent, and the dependence is significantly 
different in various limiting cases. 

We finally compare the amplitudes of the trispectra from a contact interaction diagram with those in the de Sitter 
inflation. The amplitude of the trispectra for the planar modes in all the limiting cases I, II and III can be greater 
than those in the de Sitter inflation. For the nonplanar high-momentum modes, the amplitude of the trispectrum in 
the limiting case I can be as large as that in the de Sitter inflation, while those in the other two limiting cases cannot 
be as large as those in the de Sitter inflation. Thus, combined with the direction dependence, the trispectrum in 
the limiting case I for the nonplanar high-momentum modes, as well as those in all the limiting cases for the planar 
modes, would be particularly interesting to distinguish models. 

B. Contribution from the scalar exchange interaction 

The other leading order contribution to the trispectrum is given by the scalar exchange interaction. The trispectrum 
is given by 

(0(f J ki)0(t,ka)0(t,k 8 M*,k4)) 

= - f dt 2 r^ 1 ([iif ) (< 1 ),[iif ) (i 2 ),0 / (i,k 1 )0 7 (t,k 2 ) ( /) / (t,k3)0 / (i,k4)]]), (77) 

J to J to 

where the interaction Hamiltonian is given by Eq. (jMJ) . The trispectrum in the t — > limit takes the following form 
(0(O,ki)0(O,ka)0(O,k3)0(O,k4)) = (27r)- 3 -^ r T s (k 1 ,k 2 ,k 3 ,k 4 ) ( 5(k 1 +k 2 + k 3 + k 4 ). (78) 
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After computations, we obtain 
T s (ki,k2,k 3 ,k 4 ) 



9#„ 8 



4(fcifc 2 fc 3 fc 4 ) 

Re{^i2« -A^)(A+ 2 -A^)(Ai 3 -A^YiA^-A^y 

- K 3 KAt 2 Ah, h, fcia) + A^A^A^X(-k 3 , -k 4 , -k 12 ) 
KsKA^ 1 ^ 3 ^^ - fc i2) + 2 perms) 
KAkAk^ 1 ^,-^, -k 12 ) + 2 perms) 

- (<<< 2 )*2(-fci,-fe,-fci2) + (A ki A k2 A ki2 )*I(fc 1 ,fc 2 ,fc 12 ) 

(^k^ka^ku)*^-*!! - fc 2, feia) + 2 perms) 

{KA^KS 11 ^-^ fc 2- fc i2) + 2 P erms ) 

u« - A k2 ) (A+ - A k3 ) « - A k4 )*« - A ki )* 

AiAiA^ifo, ki, -ku) + 2 perms) 
^k^k^k!! 2 ^, -fci, -fc u ) + 2 perms) 

(A+ 2 A+ 3 A k J*l(-k 2: -k 3 , ku) + 2 perms) 
{K 2 K 3 KS T (- k ^ fc 3, fcu) + 2 perms)] }. 



(79) 



Thus the trispectrum from the scalar exchange interaction depends on k\ 2 = |ki + k 2 | and fci4 ~ |ki + k 4 | as well 
as on hi (i — 1,2,3,4) because of the internal line in the diagram. Since \Z(piiP2,P3)\ is increasing in the limit of 
Pi + P2 + P3 —> 0, the bispectrum has the maximum amplitude. 

We focus on limits fci + fc 2 — > fci 2 and fc 3 + fc 4 — > fci 2 , and k\ + fc 4 — > fci 4 and fc 2 + fc 3 — > fci 4 (case I). We also 
investigate the limiting cases of k\ = fc 2 ^> fci 2 and fc 3 = fc 4 ^> fci 2 , and fci ~ fc 4 ^> fc'i 4 and fc 2 = fc 3 ^> fci 4 (case II). 



1. Case I 



In the limit of k\ + fc 2 — > fci 2 and fc 3 + fc 4 — > fci 2 , the leading order contribution to the trispectrum is given by 



T s (ki,k 2 ,k 3 ,k 4 ) 



4(fc 4 fc 2 fc 3 fc 4 ) 

x Rc{ « - A ki )(A+ 2 A k2 )(A+ 3 A^Y (A+ - A^ 

x 



Ao^k^ki 2 ^ka^ki^k^ 



Aci^kaAci: 



Ad^k 2 ^k 12 



)']} 



For the trispectrum of the nonplanar high-momentum modes, we obtain 
Re{« "^kJK ~K 3 )\ A t ~<l 



A £ 3 A ii A ki2 A k 3 A k4. A ki 2 



AnAc 2 ^k 12 



)']} 



0(^,12) 2 ^0 3 



4fc? 2 



Thus 



T s (k 4 ,k 2 ,k 3 ,k 4 ) 



16(fcifc 2 fc 3 fc 4 )fc 2 2 



(80) 



(81) 



(82) 
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The ratio to the case of the dc Sitter inflation, Eq. (|A3|) , is given by 

which is of order 41b- e?. Thus the direction dependent trispectrum can be as large as that in the de Sitter inflation 

for jj- > e" 1 . This is marginally consistent with Eq. (|55|) . 

On the other hand, the trispectrum in the same limit is given by 



Alk^hki) (1 - e- 27r « 1 ' 1 )(l - e- 27r 9i.2)(l - e- 27r ?i.3)(l - e- 27r 9i^)(l - e -2T<Zl,l2) ' 



The ratio to the case of the de Sitter inflation is given by 

M 6 kf 2 



H$ fc 6 (1 - e- 7r 9 1 - 1 )(l - e- T «i.2)(l - e - 7r «i-3)(l - e-^i.4)' 



(84) 



(85) 



'o 

which is greater than unity since M- 3> 1. This bound is consistent with Eq. ([55]) for 1 < -4/- < e" 1 . 
Similar results can be obtained for fci + k^ — > ku and k 2 + k 3 —t fc 14 . 

2. Case II 

In the limit of k\ = ki 3> k\i and ^ = fc4 > fci2, the leading order contribution from the nonplanar high momentum 
modes is given by 

T s (kx,k 2 ,k3,k4) g ° U f kl r 2 (Q(r±,i) + Q(r± l2 ))(Q(r±,3) + Q(r± A )). (86) 
16(fc?A|) 

which is of the same order as Eq. ([82]) . The ratio to the case of the dc Sitter inflation, Eq. (| A3|) , is given by 

{hk 3 )i H$r${Q(r ±ll ) + Q(r ± , 2 ))(Q(r ± , 3 ) + Q(r ±A )) 
= ^ (fclfc fc 6 lg " (Q(r±,i) + Q(^.2)) (Q(r±, 3 ) + Q(rL, 4 )), (87) 



which is of order 4tb-£*- Thus the direction dependent trispectrum can be as large as that in the de Sitter inflation for 

■jj^ > e" 1 . This is marginally consistent with Eq. ([55]) . In addition, the trispectrum in case II, Eq. ((55)) . is suppressed 

to that in case I, Eq. (|82|) . by the ratio of ^ ^ 1- Thus in contrast to the case of the bispectrum, the trispectrum 
in case I becomes greater than that in case II. 
On the other hand, for the planar modes, 



T s (k u k 2 MM" 



4(fc 2 fc|) (1 - e- 7r « 1 - 1 )(l - e~ 7r « 1 ' 2 )(l - e- 7r «' 1 . 3 )(l - e- 7r « 1 -*)(l - e'^^) 
x je-^fe.i+ei.a) cos (2^1 - 2^ 3 ) + e " n[qi - 1+qiA) cos (2^i - 2i/> 4 ) 

+ e - w ^ 2+qi ^ cos (2^2 ~ 2^3) + e _,r (9 I ' 2+9l ' 4 5 cos (2^2 ~ 2^ 4 )}- 

The ratio to the case of the de Sitter inflation is given by 

M 6 feffcf 1 

flf fc 6 (1 - e -' r< ? 1 - 1 )(l - e- 7r « 1 - 2 )(l - e- 7r ? 1 . 3 )(l - e- 7r « 1 . 4 )(l - e-^* 1 -! 2 ) 
x je-^i.i+si.s) cos (2-0i - 2^3) + e - 7r (« 1 - 1+< ? 1 ^' cos (2^1 - 2^ 4 ) 

+ e -»(«.»+«M) cos (2^ 2 - 2^3) + e _,r( ' 1 ' a+,1 '*5 cos (2^ 2 - 2V> 4 )}, (89) 
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which is greater than unity since ^> 1. This bound is consistent with Eq. (|35|) for 1 < < e„ 1 . The ratio of Eq. 



84J) to Eq. 



is given by 



e -n( qiA + qi , 3 ) pQg^^j _ 2 ^ 3 ) + e -" (<?i,i+?i,4) 008(2^! - 2^ 4 ) 



+ e - ^^ 1 - 2 ^ 1 ' 3 ) cos(2V> 2 - 2^3) + e - 7r(9l ' 2+?1 - 4) cos(2^ 2 - 2^4) 



(90) 



which can be larger than unity for k\ 2 <C k as for the nonplanar high momentum modes. The result is in contrast to 
the case of the bispectrum from the contact diagram. 

Similar results can be obtained for the case of ki = fc 4 3> fci4 and k 2 = k^^> fei4- 



3. Features of the trispectra from the scalar exchange interaction 



We summarize the features of the trispectra from the scalar exhange interaction diagram and compare them with 
those in the de Sitter inflation. In our background the amplitude of the trispectrum in the limit of case I (for example, 
for k\ + k 2 — > k\2 and k% + k± — > ki%) is greater than those in the limits of case II (for example, for k\ = k% ^> k\ 2 
and ^3 = &4 ^> ki 2 ), by a factor of where k 3> Hq is the typical size of the momentum vectors. 

We then summarize the direction dependence of the trispectra. Though we focus on the nonplanar high-momentum 
modes, the essence is the same for the planar modes. In the limiting case I, for example, for ki + k 2 — > fei2 and 
&3 + ki — > ^12), the trispectrum is proportional to <3(r^.i2) 2 and hence depends on the direction of the internal 
momentum vector k.12. In the limiting case II, for example, for ki = k 2 k\ 2 and kz — &4 3> fci2, it is proportional 
to (Q(rj_,i) + Q(r±, 2 )) (Q(r±^) + Q(rj_, 4 )) , which depends on the directions of all the momenta. Thus the trispectra 
are direction dependent, and the dependence is significantly different in various limiting cases. 

Wc finally compare the amplitudes of the trispectra from the scalar exchange interaction diagram with those in 
the de Sitter inflation. The amplitudes of the trispectra from the planar modes in all the limiting cases I and II can 
be greater than those in the de Sitter inflation. For the nonplanar high-momentum modes, the amplitudes of the 
trispectra in both the limiting cases I and II can be at least comparable to those in the de Sitter inflation. Thus, 
combined with the direction dependence, the trispectra in all the limiting cases for the nonplanar high-momentum 
modes, as well as those for the planar modes, would be interesting to distinguish models. 



VI. SUMMARY 



We have investigated the higher order spectra of a scalar field induced by the higher order time-derivative interac- 
tions, which becomes dominant above the energy scale M, in the universe with a primordial anisotropy. Along the 
lines of our recent works [ToL Hlj , a change in the vacuum causes new features in the higher order spectra. We have 
evaluated the leading order contributions to the bispectra and trispectra. 

The leading order contribution to the bispectra is given by the three-point interaction. The mixing of the negative 

frequency mode gives peaks of the bispectra in the limits of k\ — > k 2 + k^ and its permutation, where ki = |kj| are 

the magnitudes of the momentum vectors of a triangle. An imporant feature of the bispectra is that in the limit of 

&3 — > k\ + k 2 they are suppressed to that in the limit of k\ = k 2 ^> k^ by the ratio ^ , in fixing the total momentum 

k = ki + k 2 + &3. For the nonplanar high- momentum modes, for various shapes of the triangle the ratio to the case 

3 — ± 

of the de Sitter inflation is given by "jp-^i, which is larger than unity for > e* 3 , where e* is the adiabaticity 

parameter evaluated at the instance of the matching. However, this bound is not consistent with j^- < e^T 1 , which is 
obtained from the requirement that the time when the momentum of a mode becomes smaller than the energy scale 
of the derivative interactions becomes later than the time when the de Sitter mode function starts to become valid. 
For the planar modes, the bispectra are always larger than in the de Sitter inflation by a factor (^-) < e,T 3 . In 
summary, for the nonplanar high-momentum modes the bispectra cannot be as large as that of de Sitter inflation. 
On the other hand, for the planar modes, they can be greater than those in the de Sitter inflation. The bispectra are 
direction dependent, and the dependence is significantly different in various limiting cases. Thus the bispectra for the 
planar modes would be particularly interesting to distinguish models. 

We have also investigated the trispectra in the same spacetime background. The leading order contribution to 
the trispectra is obtained from two kinds of diagrams: the contact interaction and the scalar exchange interaction 
diagrams (see, e.g., fl9j for the Feynman diagrams). We have first considered the nonplanar high- momentum modes. 
The trispectrum from the contact interaction is amplified in the limit of ki k 2 + fc 3 + fc 4 (and in its permutation) , 
compared to the case of the de Sitter inflation. The ratio of the trispectrum in the limit of k% — > k\ + k 2 + k% in the 
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initially anisotropic universe to that in the isotropic one is given by j^-cf, which is larger than unity for -jj^ > e* 1 . 
This bound can be consistent with jj- < e^ 1 . Thus the trispectrum in this limit can be enhanced for smaller energy 
scales M than in the case of the bispectrum. In the other important limits, its ratio becomes of order M^-el, which 

is larger than unity for > 5 . However, this bound cannot be consistent with < e~ l . Thus it is hard to 
enhance these shapes of the trispectra. For the planar modes, the trispectra from the contact interactions are always 
larger than in the de Sitter inflation by a factor (jj^) < e^T 5 - In summary, for the nonplanar high-momentum modes, 
the trispectra can be as large as that of de Sitter inflation only for the limit of fci — > ki + ks + k^ (and also for its 
permutations). On the other hand, for the planar modes, these can be greater than those in the de Sitter inflation. 
The trispectra are direction dependent, and its dependence is significantly different in various limiting cases. Thus 
for the nonplanar high-momentum modes, the trispectra in this limiting case would be particularly interesting to 
distinguish models. Also, as for the bispectra, in all the limiting cases the planar modes would be important. 

In the case of the trispectrum from the scalar exhange interactions, the amplitude in the the limit of k\ + /c2 — > 
|ki + k:2 1 and k% + k^ — > |ki + k 2 | (and also for k\ + fc.4 — > |ki + k4 1 and k^ + k 3 — > |ki + k^) can be more important 
than that in the limit of k\ = fc 2 |ki + k2 1 and k$ = k± ^> |ki + k2 1 , by the ratio This result is in contrast to 
the case of the bispectra and would give an unique signature from an initially anisotropic universe. For the nonplanar 
high-momentum modes, the ratio of the trispectra in the initially anisotropic universe to those in the case of the de 
Sitter inflation is given by ^e*, which can be more important for > e" 1 . This bound can be marginally consistent 

with < e~ 1 . For the planar modes, the trispectra from the contact interactions are always larger than in the de 

Sitter inflation by a factor (jj^) 6 < e,T 6 - In summary, for the nonplanar high- momentum modes the trispectra can be 
as large as that in de Sitter inflation in all the limiting cases. Also, for the planar modes, they can be greater than 
those in the de Sitter inflation. The trispectra are direction dependent, and its dependence is significantly different in 
various limiting cases. Thus for the nonplanar high-momentum modes, the trispectra in all the limiting cases would 
be particularly interesting to distinguish models. Also, as for the bispectra, in all the limiting cases the planar modes 
would be important. 

Before closing this article, we would like to mention the directions of our future studies. One of the most important 
directions is to formulate the nonlinear cosmological perturbations theory in the anisotropic universe and evaluate the 
higher order spectra. Since in the anisotropic cosmological background, even at the level of the linear perturbations, 
there is a significant mixture of the different types of perturbations in terms of three-dimensional symmetry, the same 
effects may also become important during the isotropization, as well as in the change of the vacuum state. Another 
important direction is to cure the initial curvature singularity in other types of initially anisotropic spacctimc which is 
singular, and to present a successful quantization scheme of the perturbations. They may be achieved by the matter 
fields, which could significantly backreact on the spacetime geometry and regularize the initial singularity. 
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Appendix A: Higher order spectra of a scalar field in the de Sitter universe 

In this appendix, as a reference, we discuss the higher spectra of a scalar field in the de Sitter universe, induced by 
the derivative interactions given in the main text (|34l) . 



1. Bispectrum 



The bispectrum from the contact interaction in the de Sitter inflation is given by 



B(k 1 ,k 2 ,k 3 ) = - ° (Al) 

(fcifc 2 K3)(fci + h + hy 
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2. Trispectrum 

a. Contribution from the contact interaction 



The trispectrum from the contact interaction in the de Sitter inflation is given by 

79 fj s 

T c (k 1; k 2 ,k 3 ,k 4 ) = ——— (A2) 
(fcifc 2 K3Mj(fci + k 2 + k 3 + k^r 

b. Contribution from the scalar exchange interaction diagram 

The trispectrum from the scalar exchange interaction diagram in the de Sitter inflation is given by 

rp /i . , k] 9ffg f ku ku \ r .„s 

Sl U 2 ' 3 ' 4> 2{k 1 k 2 hk 4 )\(k 1 + k 2 + k 12 )3(k 3 + k 4 + k 12 )3 + (fc 2 + fc 3 + fc 14 )3(fc 1+ fc 4 + fc 14 )3/- 6 > 
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